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1. Introduction 

Discretization errors are an important source of uncertainty in lattice QCD calculations with 
heavy quarks. As with gluons and light quarks, one uses effective field theory to understand and 
control discretization (aka cutoff or lattice-spacing) effects. With heavy quarks there are two pos- 
sibilities: a (modified) Symanzik effective Lagrangian and a (modified) heavy-quark effective La- 
grangian (HQET for heavy-light hadrons; NRQCD for quarkonium). 

Usually, the Symanzik local effective Lagrangian (LEJzf) is devised with only one short dis- 
tance in mind, the lattice spacing a. Similarly, HQET/NRQCD is usually constructed with only the 
heavy quark's (or quarks') Compton wavelength(s), 1 jrriQ, as short distance(s). With heavy quarks 
on a lattice one must consider both kinds of short distances. In particular, couplings in the lat- 
tice Lagrangian and Wilson coefficients in the effective Lagrangians depend on the dimensionless 
ratio(s) m,Qa. 

Treating the inverse quark mass as a short distance, the Symanzik LEJ2? can be written ^ 

^Sym = -q\ Y4D4 + J—y-D + mAq + aK] al q{y A D 4r + m Y ) 2 q + aKf t qD 2 q 
\ V m 2 J 

+ aK^qiL ■ Bq + aK^qCC -Eq+---, (1.1) 

where m\ is the renormalized rest mass, and mass-dependent short-distance effects are lumped into 
\Jm\jm2 and the coefficients K M . An improved action is devised by adjusting Jz?sym to reproduce 
QCD. Some coefficients may be made to vanish via field redefinitions, such as and K\ at [§ g]. 
The others must be addressed by matching the underlying lattice Lagrangian: 

mi =m 2 = m Q , 4 at (c B )=0, 4 at (c £ )=0,..., (1.2) 

where iuq is the physical quark mass, and cb and ce are couplings of dimension-five interactions 
in the Fermilab lattice action for heavy quarks [^]. 

Treating the lattice spacing and inverse heavy quark mass as short distances, the heavy-quark 
effective Lagrangian is [Q] 

-2hq = -h {±) (74O4 + mi) h& + — + S- — 

where the Z lat capture the short-distance effects. In parallel, one can construct this effective La- 
grangian for continuum QCD, with the same operators but different coefficients Z cont . Now im- 
provement is achieved by choosing improvement couplings in the underlying lattice Lagrangian, 
so that the two effective HQ Lagrangians coincide. That is 

m 2 = m Q , Zf{c a ) = Z- n \ Zf{c E ) = Z™ nt , • • • • (1.4) 



Both sets of conditions, Eqs. (1.2) or Eqs. (1.4), yield the same results for the improvement cou- 
plings cb and c#. 
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In this paper we focus on an extension of the Fermilab method to include all interactions of 
dimension six. Some aspects have been reported earlier []|, gj. Here, we first discuss a subtlety 
in the power counting, which implies that we must also consider certain interactions of dimension 
seven to be commensurate with some of those of dimension six. 



2. Power Counting 

With heavy quarks one must pay special attention to power counting. This is simplest in 
heavy-light hadrons, where the typical three-momentum in the rest frame is p rs Aqcd- Then it is 
useful to define nr, which is or 1, depending on whether the Dirac matrix T in a quark bilinear 
commutes or anticommutes with y 4 . A bilinear of dimension d introduces physical effects of order 

which are all the same once one allows for coefficient functions that depend on raga, with niQa 
of order unity. This power counting has recently been considered by Christ, Li, and Lin It is 
useful for hadrons with one heavy quark, but for quarkonium one should adopt the power counting 
of nonrelativistic QCD ||. 

We shall denote terms in the lattice action by „ r ) to classify them by their power count- 
ing (2T ). The terms with « r = 1 are necessary to ensure a smooth limit when niga — > [|[], which is 
a feature distinguishing the Fermilab method from lattice NRQCD. Ref. ^ treated 5( 5 ) and 5(5,1) 
(and cursorily S^fl))- Here we complete the analysis presented at earlier Lattice symposia [j|, ||] to 
encompass the full set of interactions in 5(6,1) and 5(7,o)- 

3. Improved Fermilab Action 

The Fermilab action is a generalization of the Sheikholeslami-Wohlert action. We write 

l 

5 = 5 +^ J! *V,nr)> ( 3 -l) 
d=5 «r=0 

where 

S = m()a 4 ^^)^(x)+a 4 ^^(x)74D 41at VA(^)-^a 5 ^^(x)A4 lat VA(x) 

XX X 

+ C« 4 £ V/(x)y-D lat ^(x) - ir^L^MA^VW. (3.2) 

X X 

We denote lattice fermions fields with y to distinguish them from the continuum quark fields in 
Eqs. (1.1) and (1.3). The coupling £ allows one to adjust ni2 = m\ in Eq. ( |1.1[ ). The notation for 



Di a t, A lat , etc., is as in Ref. [^J. The dimension-five interactions are [§, ||] 

5(5,o) = S B = -\c B £a 5 Y,Y(x)iZ-BMY{x), (3.3) 

X 

5(5,1) =S E = -\c E ^a 5 Y,W) a - E ^M^)^ (3.4) 
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where the notation Sb and S E is from Ref. [|3j]. At dimension six and seven we introduce 
%o) = r E a 6 J^y(x){Y ■ D M ,a ■ E M }y(x) 

X 

+ Z E a 6 Y,Y(x)74 (Olat ■•Blal-^lat -Oiat) (3.5) 

x 

5(6,1) = Cl^LWL^AlatAiiatV^W 

+ ^^^Mr-AatA^M*) 
+ c 3 a 6 £i//(^){y-Di at ,/E-Bi at }i//(x) 

+ Z 3 a 6 £ V/(x)y • (Aat X Slat + Si at X Aat) Y(x) 

x 

+ c EE a 6 Y,¥i.x){r4D4uuOC ■ E m }y(x), (3.6) 

x 

5(7,o) = c^^v/W^A^V^W 

x i 

+ c 5 a 7 £ \fr(x) ££{^131, A hat }y(x) 

+ r 5 a 7 £ y(x) £ £ /£,• [A-BAjU V (x) 
x ' 

+ z 6 a 7 £^(x)(A, ( a 3 t ) )Vw 

+ zya^^lAj^/E-Biatj^x) 

x 

+ 4fl 7 £^(x)[Dia:-aA]iatvW 

X 

+ r 7 a 7 £ y/(x) y • Aat • Bi at y • Aat y(x) 

x 

+ r'y £ i//(x) [D • (fi x D)]iatVA(x) 

x 

+ r BB a 7 £ yr(x) (zX • Bi at ) 2 yf(x) 

X 

+ ZBB« 7 £^WBlafBiat^(x) 

x 

- £ VK*) (« • Slat) 2 W(x) 

x 

+ z ££ a 7 £ \j/(x)E u ■ E u y(x)- (3.7) 

X 

In particular, in heavy-quark power counting the interactions in S^o) are commensurate with those 
in 5(6 i). Most of these interactions have been discussed previously [|5[ ||], but the following is the 
first discussion of interactions with two chromomagnetic or chromoelectric fields, and the related 
interaction with coupling c EE . 1 

By applying arbitrary field redefinitions one can make certain coefficients in the effective La- 
grangian vanish, and their interactions are called redundant. A corresponding number of couplings 

'At dimension six and higher one must also consider four-quark interactions, which we shall discuss elsewhere. 
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in the lattice action may be chosen by convenience. For example, to obviate the doubling prob- 
lem one sets r s ^ and r t = 1 in Eq. (3.2). One can show that all interactions with further time 



derivatives acting on quark and antiquark fields are redundant, so we omit them from Eqs. ( |3.3| )- 
([O]). Among the rest there is some freedom on which ones to call redundant; with an eye towards 
efficiency in computing quark propagators we choose 

r E = r 5 = r 1 =r' 1 = r BB = r EE = 0, (3.8) 

a choice that is substantiated further by the matching calculations discussed in the next section. 



4. Matching 



To find out how to adjust the other couplings, we must explicitly derive Eqs. (1.2) or, equiva- 



lently, Eqs. (1.4), and solve them for the couplings in the lattice action. Because the action includes 
interactions with two chromoelectric and two chromomagnetic fields, we have done so by working 
out the amplitude for Compton scattering. We find, not unexpectedly, that many couplings vanish 
at the tree level. To be specific, tree-level matching yields 



ZE = Z3= Z(, 



Zl 



z'q = z,bb = z EE = 0, 



(4.1) 



but these become non- trivial beyond the tree level. 

The other couplings — cb, ce, cee, c\, c%, C3, C4, and C5 — must be non-zero, with an explicit, 
sometimes nontrivial, dependence on (moa, £,r s ). We do not have space to discuss the intermediate 
steps here, so we quote the final result of tree-level matching. We find 



cb 
ce 

c\ 

C2=C 3 



C 2 -i 



< + ■ 



r s £ rjmoa(2+moa) 



moa(2 + moa) ' \+moa 
moa(2 + moa 



+ 



4(1 +m a) 2 



1 y ' 

6(l+mofl) 



(4.2) 
(4.3) 

(4.4) 



C 3 (C 2 -i) 



C 2 [C + 2r s (l +m a) -3r. v £ 2 /(l +m a)] 



[2nioa(2 + niQa)Y 



8moa(2 + mna) 



+ 



c E E [2 + m a(2 + mod)] 



16(1- 



■ + 



moa ) 



moa(2 + rao<2)r 2 £ 
32(1 +m a) 2 



1 + niQa 



1 



(4.5) 



l)(l+moa) c E ^ 2 +m Q a) 



+ 



[moa(2 + moa)] 2 
C(r s £-l-moa) 



+ 



C4 
C5 



2moa{2 + moa) 



+ 



mna(2 + moa) 
\r s CE^ 2 - jclC(l +m a), 



(4.6) 

(4.7) 
(4.8) 



Some aspects are not surprising: for example C2 = C3, which makes sense because together they 
provide an interaction of the form \j/(y ■ D) 3 y. Two out of the three couplings (04,05,^) must be 
nonzero, essentially to correct errors from the Bi at in 5(5,0) • We choose C4 7^ because we expect 
its interaction to be easier to compute than that with coupling . 
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5. Outlook 



In several recent calculations with Fermilab quarks, the largest systematic uncertainty comes 
from a (conservative) estimate of heavy-quark discretization effects. With the improved action 
presented here, the same technique for estimating the uncertainties ^ suggests that these effects 
are reduced to a few percent To achieve this target, tree-level matching should suffice for S(6,i) 
and 5(7 0). F° r ^(5,0)» ^(5,i)> anc ^ ^(6,0) higher accuracy is needed, either at the one-loop [10, 11] or 
(for cb in S^ o)) nonperturbative level fll2|]. 
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